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We consider quasi-stationary scattering of interacting bosonic matter waves in one-dimensional 
waveguides, as they arise in guided atom lasers. We show how the truncated Wigner (tW) method, 
which corresponds to the semiclassical description of the bosonic many-body system on the level of 
the diagonal approximation, can be utilized in order to describe such many-body bosonic scattering 
processes. Special emphasis is put on the discretization of space at the exact quantum level, in 
order to properly implement the semiclassical approximation and the tW method, as well as on the 
discussion of the results to be obtained in the continuous limit. 


I. INTRODUCTION 

The perspective to realize atomtronic devices M 
as well as the exploration of transport features that 
are known from electronic mesoscopic systems a a 
have strongly stimulated the research on the dynamical 
properties of ultracold atoms in open systems. While 
fermionic atoms provide direct analogies with the elec¬ 
tronic case he the use of a bosonic atomic species 
brings along new aspects and challenges for the atomic 
transport problem ||, f|, related, in particular, with the 
link between a mesoscopic Bose-Einstein condensate in a 
reservoir and the microscopic dynamics of an ensemble of 
few interacting atoms within a transistor-like device. A 
particularly promising configuration for the experimen¬ 
tal study of these aspects is provided by the guided atom 
laser 10l4l2j in which atoms are coherently outcoupled 
from a trapped Bose-Einstein condensate into an opti¬ 
cal waveguide. A coherent atomic beam can thereby be 
created and injected onto engineered optical scattering 
geometries. This would allow one to study bosonic many- 
body scattering at well-defined incident energy. 

A theoretical modeling of such waveguide scattering 
processes within guided atom lasers faces the challenge 
of dealing with an open system in a many-body context, 
potentially involving a very large number of atoms in 
total. Exact diagonalization methods therefore quickly 
encounter limitations when describing such processes. 
The Gross-Pitaevskii approximation, on the other hand, 
which has been used in a number of studies on guided 
atom laser processes ^3j-[l7j. is granted to work in the 
mean-field limit of large atom densities and vanishing 
atom-atom interaction strengths [l8j]. It is, however, 
known to break down in the presence of significantly 
strong interaction effects 0] . 

The truncated Wigner (tW) method f20l - [23) appears 
to be a reasonable compromise in this context. Accord¬ 
ing to common understanding, this method essentially 
consists in representing the time evolution of a bosonic 
many-body quantum state by classical fields evolving ac¬ 
cording to a Gross-Pitaevskii equation, whose initial val¬ 
ues are chosen such that they properly sample the ini¬ 
tial quantum state under consideration. They thereby 


provide a stochastic sampling of the approximated time 
evolution of the many-body Wigner function defined in 
the phase space of the bosonic field. The tW method has 
been successfully applied in a number of contexts involv¬ 
ing dynamical processes of ultracold bosonic atoms. This 
includes, most recently, scattering of atom laser beams 
within one-dimensional Bose-Hubbard systems (24}. In 
this latter transport context, the tW method was shown 
to yield quantitative predictions for the average trans¬ 
mitted current that were shown to be in good agreement 
with complementary matrix-product state calculations in 
the regime of moderate on-site interaction strenghts [24[ . 

While it is possible to formulate the tW method in 
a functional manner [25), its practical implementation is 
most conveniently achieved by means of a discrete basis of 
the single-particle Hilbert space. This is straightforward 
to accomplish for dynamical process that are effectively 
taking place within spatially confined regions. In that 
case, periodic boundary conditions can, e.g., be imposed 
at a sufficiently large distance from the quantum many- 
body wave packet to be studied, giving thereby rise to an 
effective discretization in momentum space. This is, how¬ 
ever, not a viable option for studying quasi-stationary 
scattering processes which are generally characterized by 
an infinite spatial extension. We therefore propose to 
discretize the position space in order to implement the 
tW method in this context, as was effectively done in 
Ref. [24) A primary purpose of this paper is to introduce 
this discretization procedure in some detail and discuss 
its validity in the continuous limit, both from an analyt¬ 
ical and from a numerical point of view. 

Moreover, we present in this paper an unconventional 
derivation of the tW method in the framework of the 
semiclassical van Vleck-Gutzwiller theory |26j. This 
allows us to identify the tW method with the diag¬ 
onal semiclassical approximation in the bosonic field- 
theoretical context, and, as shown in Ref. [27), to quan¬ 
titatively account for interference effects beyond tW. We 
note that the tW method can also be derived through 
other approaches such as Wigner-Moyal expansions |21- 
f23[ | . quasiclasical corrections to the effective action |28) , 
and the so-called semiclassical approximation in the con¬ 
text of the Keldysh approach [29). In those approaches, 
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however, quantum corrections are perturbatively incor¬ 
porated on top of a classical background given by the 
propagation of phase space distributions along classical 
trajectories that are accounted for in an independent (i.e. 
incoherent) manner. Interference effects involving differ¬ 
ent trajectories are essentially non-perturbative and re¬ 
quire special resummation techniques that are justified 
only in the presence of small parameters, (typically the 
strength of interactions or the size of quantum fluctu¬ 
ations). While the van Vleck-Gutzwiller theory yields 
exactly the same (tW) approximation to leading (classi¬ 
cal) order as the other approaches mentioned above, its 
key asset resides in the fact that it can readily incorpo¬ 
rate quantum interference in a non-perturbative manner 
and describe, e.g., coherent backscattering in the Fock 
space of many-body systems [23]. This is our primary 
long-term motivation to advertise this approach in this 
article. 

Our paper is organized as follows: In Section[TT]we shall 
discuss in some detail the spatial discretization procedure 
of one-dimensional waveguide scattering configurations. 
Section m is devoted to deriving the tW method in the 
framework of the van Vleck-Gutzwiller theory. We shall 
then argue in Section IIV Al that the tW approach can 
be reformulated in terms of a stochastic Gross-Pitaevskii 
equation whose noisy components exhibit well-defined 
statistical characteristics in the continuous limit. These 
findings are confirmed by numerical results on transport 
through a symmetric double barrier potential, as we show 
in Section IIVBI Both the average total current of atoms 
and its incoherent part tend to finite values in the con¬ 
tinuous limit. 

II. DISCRETIZATION OF SPACE 

We consider a many-body scattering process of a coher¬ 
ent atomic matter-wave beam within a waveguide. This 
matter-wave beam is supposed to be created by a co¬ 
herent outcoupling process from a trapped Bose-Einstein 
condensate, as is commonly done in guided atom lasers 
EMI- If we assume that only the transverse ground 
mode of the waveguide is populated, we can describe this 
system by the many-body Hamiltonian 

r°° „ / t2 n2 'v 

H = — — + 

l - 

+— / g(x)ip^ (x)^ (x)ip(x)ijj(x)dx 

+ f k(x) ($(x)b + &'<p(x)\ dx + fjPb (1) 

J —OO ' ' 

where i/j'(x) and ^j(x) respectively represent the creation 
and annihilation operators of a bosonic particle at the 
longitudinal position x within the waveguide, m is the 
mass of the atoms, V(x) describes a scattering potential 
(given, e.g., by a sequence of barriers) within the waveg¬ 
uide, and g(x) is the spatially dependent one-dimensional 


interaction strength of the atoms in the wave guid e. It is 
approximately given by g{x) ~ 2Hui±(x)a s (x) [301 ] where 
w_l is the transverse confinement frequency of the waveg¬ 
uide and a s is the (generally very small) s-wave scattering 
length of the atomic species under consideration. Both 
co± and a s can depend on the position in the waveguide, 
the former through a spatial variation of the transverse 
confinement, and the latter through a spatially depen¬ 
dent Feshbach tuning. 

The trap is modeled by a single one-particle state with 
energy /i whose associated creation an annihilation opera¬ 
tors are given by ¥ and b. Trapped atoms are outcoupled 
to the waveguide through the spatially dependent (real) 
coupling strength n(x), which, e.g., in Ref. would 
model the effect of a radiofrequency field that flips the 
spin of the atoms. We shall in the following assume a 
macroscopically large initial population N —> oo of the 
atoms in the trap in combination with a very small out- 
coupling strength k —> 0, such that the trapped conden¬ 
sate is not appreciably affected by the outcoupling pro¬ 
cess on finite evolution times. The time evolution of the 
system can then be effectively described by the equation 

+V(x) - pj i)(x,t) + K(x)be~^ t/h 
+g(x)ft(x,t)ip(x,t)'ip(x,t) (2) 

for the time-dependent field operator ip(x,t) describing 
atoms in the waveguide. 

In order to implement the tW method for this many- 
body scattering problem, we first need to introduce a dis¬ 
cretization procedure of this spatially continuous quan¬ 
tum field equation. As pointed out above, such a dis¬ 
cretization cannot be defined in momentum or energy 
space (e.g. through the introduction of periodic bound¬ 
ary conditions at x = ±L for some large L — > oo) as 
this would not be compatible with the formation of a 
quasi-stationary scattering state. We therefore propose 
to discretize the position space, namely through the in¬ 
troduction of a high-energy cutoff in momentum space 
at p = ±7 vh/S for some effective spatial grid size 6 —> 0. 
Correspondingly, we also modify the Hamiltonian of the 
free kinetic energy close to the high-energy cutoff, such 
that it reads 

/ OO 

$(x) (1 — cos (Sp/h)) r ij){x)dx (3) 

-oo 

with Eg = h 2 /(mS 2 ) and p = —ihd/dx. This Hamil¬ 
tonian is still diagonalized in the eigenbasis of the 
normalized waves 4>k{x) = exp(ifcx)/\/27r satisfying 

f_Z ^k( x )^k'(x)dx = S(k — k'). Their associated eigen¬ 
values now read Ek = Eg(l — cos kS) close to the cut¬ 
off while they are still approximately given by Ek — 
h 2 k 2 /(2m) for kS <C 1. 

In close analogy with the theory of spatially periodic 
systems, we now introduce an effective Wannier basis 
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through the spatially localized functions 


M x ) 


JIf 

V 2-7T J_ 

'7r/<5 

(j)k{x)e~ llkS dk 

-n/S 

(4) 

1 

/ 7T . 

(5) 

^=sinc 



for integer l € Z, with sinc(x) = sin(a;)/a;, which can ef¬ 
fectively be seen as Fourier series coefficients of the waves 
ifik- They therefore satisfy f^° <j>*(x)(j>i>(x)dx = 5u> and 
thereby form an orthogonal basis set that spans the re¬ 
stricted Hilbert space obtained after the introduction of 
the momentum cutoff. Defining the corresponding cre¬ 
ation and annihilation operators m , ifi such that we have 

i/j(x) = YmZ-oo i we can now reformulate the free 

kinetic Hamiltonian 0 such that it reads 


^ uu 

Ho = E (sV’/V’z - i’l+ii’i - 1) ■ ( 6 ) 

l =—oo 


This expression is identical to the one that would be ob¬ 
tained from a finite-difference approximation of the ki¬ 
netic energy in the waveguide. 

Altogether, we thereby obtain in the limit S —> 0 the 
effective Bose-Hubbard Hamiltonian 


H 


= E 




l =—OO L 


(i’i+i'4’1 + 1 ) 


Kl 


b + b^ipi 


(7) 


where we introduce the definitions V) = V(IS ), gi = 
g(l5)/5, and Ki = V5k(15), and where, for the sake of 
convenience, we redefine in Eq. 0 the zero of the energy 
scale such that it coincides with the energy of trapped 
atoms. The time evolution of the discrete field operator 
ipi is then given by 

= (E s - 11 + Vi)ipi{t) - E 

+gi${t)'4>i(t)'ipi{t) + Kib. (8) 


III. SEMICLASSICAL DERIVATION OF THE 
TRUNCATED WIGNER METHOD 

Having accomplished the discretization of the exact 
quantum description, a semiclassical approach can be 
used in the regime of large particle numbers N 1, 
which does not resort to the numerical solution of the 
operator-valued equations (El). This approach has three 
levels of approximation. The first level is the purely 
classical limit consisting of the propagation of the (dis¬ 
cretized) Gross-Pitaevskii equation that is obtained from 
Eq. ® by the substitution f)i(f) A ipi(t) and A 

if* ( t ) where the latter classical fields are canonically con¬ 
jugated variables. In a second step, a finite initial phase- 
space distribution of the quantum many-body Wigner 


function is accounted for. This phase-space distribution 
can contain non-classical correlations which is, however, 
not the case in our present problem. Finally, a full- 
fledged semiclassical approximation is obtained by co¬ 
herently adding amplitudes associated with classical tra¬ 
jectories, in order to address quantum interference ef¬ 
fects [3l|. The second stage in this hierarchy of levels 
of approximation is the tW method. Here we derive this 
method in an unconventional manner which is based on a 
semiclassical approximation of the Feynman propagator 
of our effective Bose-Hubbard system 0 in terms of a 
coherent sum over solutions of the mean-field equations. 

The notion of the term ” semiclassical” in our approach 
is conceptually different from its meaning in the context 
of the Moyal-Wigner expansion used in quantum optics 
M, the quantum corrections to the effective action 
of Ref. {2f|, and the quasiclassical corrections within the 
Keldysh approach [231. I n those approaches, the clas¬ 
sical limit is identified in terms of the transport equa¬ 
tions of probability distributions in phase space (the so- 
called classical field in the Keldysh language [29l|i. while 
quantum corrections are systematically incorporated as 
perturbation series around this classical background. Al¬ 
though interference effects, with their charateristic non- 
perturbative dependence on H , can be obtained within 
those frameworks by special resummation techniques, 
these techniques are in principle not suitable to tackle 
classically chaotic situations where all relevant physical 
forces governing the dynamics of the system under con¬ 
sideration are of the same order and no small parameter 
can be identified. 

The semiclassical van Vleck-Gutzwiller approach, on 
the other hand, is particularly powerful precisely in this 
chaotic regime where it allows one to predict universal 
interference effects. Its key feature is that it approx¬ 
imates the quantum time evolution of the system by 
a semiclassical superposition of amplitudes and there¬ 
fore explicitly accounts for quantum mechanical interfer¬ 
ence through coherent double sums over classical paths. 
By working with amplitudes instead of probabilities, we 
can explicitely identify tW with the assumption that 
pairs of different classical solutions have uncorrelated ac¬ 
tions. Going beyond this so-called diagonal approxima¬ 
tion and taking into account systematic off-diagonal con¬ 
tributions, the semiclassical approach has been used to 
predict interference effects beyond both tW and its qua¬ 
siclassical corrections in good agreeement with numerical 
simulations (27ll3^. 

Following Ref. {23, we begin with the time evo¬ 
lution of an arbitrary many body state, |(i)) = 

C/(t)|\E , (0)), governed by the time evolution operator 
U(t) = exp {—itH/h) associated with the Hamiltonian 
0. Our Hilbert space 

n = span{|n) = |... ,n_i,n 0 ,ni,.. .),m e N 0 } (9) 

is spanned by the Fock states associated with the site ba¬ 
sis, where the discretized field operators act in the usual 
manner, namely by raising or lowering the number of 
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particles in a given site: 

i>j\■ ■ ■ ,ni,...) = y/ni + 1| ... ,m + 1,...), (10) 

= y/ni\... ,m - 1,...). (11) 

The starting point of our approach is to use of a dif¬ 
ferent basis for R, which is the formal equivalent of the 
position representation used in the usual derivation of 
the Feynman propagator in first-quantized systems. This 
new basis |q) is given by the common eigenstates of the 
commuting set of hermitian operators qi = 
which are known in quantum optics as quadratures. Us¬ 
ing the defining property <?;|q) = qi\q), the orthonor¬ 
mality and completness relations (q|q') = <S(q — q') and 
/ dq|q)(q| = 1 of the quadrature basis follow directly. 

Having at hand a complete basis with continuous 
states, a path integral representation of the time evo¬ 
lution operator can be obtained in the usual way. Con¬ 
trary to the usual kinetic-plus-potential Hamiltonian in 
the first-quantized case, the use of the eigenstates |p) 
of the momentum quadratures pi = —i^/h/2(ipi — ■ipj) 
is not only useful but essential to get the amplitude 
A'(q f ,q‘,t) = (q f |U(t)|q 1 ) of the propagation from an 
initial to a final quadrature state 

K (q f , q 1 , t) = J 2?[q(s),p(s)]e lfl[q(s) ’ p(s)1/fi (12) 


as a sum over all solutions of the classical shooting prob¬ 
lem (up to an extra phase z/x 7 7t/ 4 involving the Maslov 
index /z 7 , which is not important in the following). 

The time evolution of the expectation value (A) t = 
('F(f)|A(q,p)l'F(i)} of a properly (Weyl) ordered opera¬ 
tor A = A(q, p) is given by 


(A)t = J dq f dq i dq' i (q' i |^(0))(tk(0)|q i ) 

X K (q 1 , q 1 , t)*A( q f , -ihd/dq 1 )K (q 1 , q'\ i)(16) 


Our goal is to obtain a semiclassical approximation for 
this expression. As first step, we substitute the semi¬ 
classical propagator m in Eq. m- At this level of 
approximation, we only assume that a typical action in 
Eq. (1131) is significantly larger than h. In that limit, for 
any smooth functions F(x),G(x), H(x) we have 

F (~m£j G( q )e^ = F G( q ) e ^+0(S). 

q ■ q . 

Using dR 7 ( q f , q' 1 , f)/dq f = p 7 (q f ,q‘,t), we obtain then 


A(q‘, -ihd/dq^K ^, qi) 


* E^p 


det 


d 2 i? 7 (q f ,q i ,t) 


dq‘<9q f 


2nh 


1/2 


JR -,(q ,q‘ 


in the form of an integral over paths q(s), p(s) satisfying 
the boundary (shooting) conditions q(s = 0) = q 1 and 
q(s = t) = q 1 . Here, the action functional is given in its 
Hamiltonian form, 


%(«),P(s)] = [ ds[q{s)-p(s)-H c \q(s),p{s))} (13) 
Jo 

where the classical Hamiltonian H cl is obtained from the 
quantum one by the substitutions H > ipi and H > ip* 
with ipi = (qi + ipi)/y/2h. Properly taking care of the 
Weyl ordering of operators, which yields the replacement 
ipl'ipj + 4’J'ipi H > 2 ipfipi, we then obtain from Eq. dTJ) the 
classical Gross-Pitaevskii-type Hamiltonian 


H c\ 


E 


{Eg - p + Vi)ipfipi 


Y {i’t+ii’i + i’ti’i+i) 


+ -2 ) + Ki (ipfb+Fipi) 


(14) 


up to constant terms that are not important. 

We now continue with a stationary phase analysis of 
the path integral, following Gutzwillcr’s pionering work 
in the 60’s and 70’s [26|. This programme was accom¬ 
plished in Ref. {27} • It yields as final result the propaga¬ 
tor 




det 


d 2 R 7 (q f ,q‘,t) 


dq'dq 1 


2nh 


1/2 


a j-R-,(q f ,q\t)/ft 


(15) 


This then yields a double sum over classical trajectories 
in Eq. (1161) . It is remarkably accurate even for moder¬ 
ate values of the classical actions and for arbitrary large 
times, and it sucessfully describes interference phenom¬ 
ena through the coherent sum over the oscillatory ampli¬ 
tudes that are associated with each trajectory. 

In a further approximation, one makes use of the fact 
that the double sum over paths contains a large number 
of terms that effectively cancel each other in the presence 
of an average over initial and final positions as in Eq. (1161) . 
Only pairs of trajectories with similar actions yield non¬ 
vanishing contributions due to phase cancellations, which 
gives rise to incoherent sums of slowly oscillatory terms 
with non-zero average. This is the essence of the diago¬ 
nal approximation, the standard tool to obtain leading, 
classical-like, contributions from coherent double sums 
representing semiclassical amplitudes. In the present 
context, the diagonal approximation is implemented by 
pairing 7 (q f ,q\f) and 7(q f ,q / \t) for which q 1 ~ q' 1 . In 
the semiclassical limit, all smooth terms in the propa¬ 
gator are then evaluated at q = (q 1 -(- q n )/2, while the 
actions are expanded up to first order in q 1 q' 1 according 
to R 1 {q i ,q i ,t) -i? 7 (q f ,q'\t) ~ -ply(q f ,q, t) ■ (q^q' 1 ). 

The double sums are then transformed into a single 
sum yielding the diagonal approximation for (A) t as 


(A)t = /dq f dqW i (qi^(0))(^(0)|q i )^A[q f ,pf(q f ,q,t)] 

(19) 


det d 2 Rry(q { , q, t) 


dqdq 1 


2nh 


e ->p)(q ■q.i)'(q I -q")/s 
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This expression can be further simplified by noticing that 
9 2 i? 7 (q 1 ,q, t)/dqdc^ = — dp 7 (q f ,q, t)/dq l is the Jaco¬ 
bian matrix of the transformation q f —> p 1 , as these vari¬ 
ables are functionally dependent through the solutions 
q f = q‘(q, p 1 , £), p f = p f (q, p‘,i) of the classical equa¬ 
tions of motion. We can then change the integration 
over final coordinates by an integration over initial mo¬ 
menta. Together with the substitution (qbq' 1 ) i->- (q, x) 
with q = (q 1 + q' 1 )/2 and x = q 1 — q' 1 (and with the 
replacement p 1 = p), this finally yields 

(^)t = / dqdpA[q f {q,p,t),p f (q,p,t)]Wy(p){q,p), 

( 2 °) 

where we introduce the Wigner function associated with 
the state |$) through 

w*{ q,p) = (n/S) e- lp x/ft (q+x/2|$)(<l>|q-x/2) 

( 21 ) 

with the product going over all single-particle sites 
!33l . It is possible to reformulate Eq. as 

= J dqdpA(q, p)W* (0) (q, p, t) (22) 

where W^( 0 )(q, p,t) is the time-dependent Wigner func¬ 
tion which evolves according to the well-known trun¬ 
cated Wigner equation. This demonstrates that the tW 
method is obtained by neglecting off-diagonal terms in 
the semiclassical calculation of time-dependent expecta¬ 
tion values of quantum observables. 

For the waveguide scattering configuration under con¬ 
sideration, the classical field equations to be numerically 
integrated in the framework of the tW method are de¬ 
rived from Eq. m and read 

= ( E s — M + Vi)ijji(t) - ^ (^j-i(t) + 4>i+i(t)) 

+gi (|'*M i )| 2 - l) + VNki . (23) 

Here we represent, in the limit of a large number of atoms 
TV —^ oo, the Bose-Einstein condensate in the trap by a 
coherent state with the amplitude y/N and with van¬ 
ishing relative uncertainty in its amplitude and phase. 
The initial quantum state in the waveguide, on the other 
hand, is a vacuum state that has to be correctly sampled 
in the framework of the truncated Wigner method. More 
specifically, the initial field amplitudes in the discretized 
waveguide have to be chosen as ipi(0) = (Ai + iB {)/2 
where Ai and Bi are independent Gaussian random vari¬ 
ables with zero mean and unit variance. As a conse¬ 
quence, each site l of the grid representing the empty 
waveguide has initially the average classical density 

hMOF = 1/2 (24) 

where the overline denotes the tW average over the ran¬ 
dom variables. This pseudo-density of half a particle per 


grid point is to be subtracted when calculating the av¬ 
erage atomic density p(x,t) = {^(x, t)^(x, t)) using the 
tW method: this average density is evaluated as 

p{x = IS, t) = i (l^(i)l 2 - 1/2) ■ (25) 

The current j{x,t) = Re[(i p^(x,t)p'i/j(x,t))]/m, on the 
other hand, is not articifially affected by the vacuum fluc¬ 
tuations; it is evaluated as 


j(x = 18, t) = - V’f+i ■ 

( 26 ) 

For both the density and the current we can define 
coherent components according to 


^coh / 


r h (is,t) = 


Mt) 

h 


(27) 


2 imS 2 


(V’fW ^i+i(t) - iPt +1 (t) iMi)) 2 ^) 


which result from a direct tW average of the complex 
amplitudes i/ji ■ They would correspond to the mean-field 
predictions obtained by the classical Gross-Pitaevskii ap¬ 
proach in the limit of large densities and vanishing atom- 
atom interaction strength. The associated incoherent 
components are defined through 

p incoh (l5,t) = p(l8,t)-p coh (l5,t), (29) 

j incoh (l8,t) = j{l5,t) — j coh (lS,t). (30) 


They approximately characterize the amount of quantum 
depletion in the waveguide. 


IV. THE CONTINUOUS LIMIT 
A. Quantum noise 

The tW method is certified to be valid in the semi¬ 
classical regime where the average population per single¬ 
particle state is large compared to unity. However, this 
condition is no longer met in the continuous limit of 
vanishing grid spacing S —> 0. Indeed, assuming that 
the scattering process under consideration gives rise to 
a well-defined average atom density p{x, t) in the waveg¬ 
uide, the corresponding discrete mean- field am plitudes 
would essentially be given by ipi(t) — y/Sp(lS,t). In the 
limit S —» 0 they would thereby become negligibly small 
as compared to the vacuum fluctuations which still fea¬ 
ture half a (pseudo-)particle per grid point according to 
Eq. HMD. 

The purpose of this section is not to discuss the validity 
of the tW method in the continuous limit (which would 
obviously require the comparison of numerical results 
with a complementary, and genuinely quantum, method, 
such as exact diagonalization) but to argue that it yields 
consistent results in this limit. We shall, for this pur¬ 
pose, introduce a Bogoliubov decomposition of the tW 
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amplitude according to ipi (t) = xi (t) + Vi (t) where \i re¬ 
sults from the unperturbed time evolution of the vacuum 
fluctuation according to 

ih §t Xl ^ = ^ Es ~ M + Vl ^ Xl ^ ~ ir te-iW + Xi+i(t)) 

+9i (\xi{t)\ 2 - 1) Xi(t) (31) 

with the initial condition Xz(0) = ^>;(0). This equation 
will not change the characteristics of the vacuum fluctu¬ 
ations in the course of time evolution as it obviously cor¬ 
responds to the truncated Wigner modeling of an empty 
waveguide. Through Eqs. (l23l) and m, Vl = ipl- Xi 
then evolves according to 

= (Eg - n + Vi)vi(t) - ^ 

+9i\Vl(t)\ 2 Vl{t) + VNki 

+9i [(2|Xi(*)| 2 - +xf(t)<PKt)] 

+9i [(2Xi{t)\<Pi{t)\ 2 +X*{t)<Pi(tj] . ( 32 ) 

which effectively corresponds to a stochastic Gross-Pita- 
evskii equation. Indeed, the first two lines of Eq. (l32l) 
describe the mean-field Gross-Pitaveskii modeling of the 
waveguide scattering process, with <^(0) = 0 as initial 
condition in the waveguide. The other two lines contain 
quantum noise terms that randomly perturb the evolu¬ 
tion of the coherent field ipi emitted by the outcoupling 
process from the trap. 

It is instructive to quantitatively evaluate the charac¬ 
teristics of this quantum noise in the continuous limit 
5 —> 0. In that limit, we can neglect the last line of 
Eq. (152)1 since this particular noise term is suppressed to 
the noise terms in the third line of Eq. ED by a scaling 
factor VS (noting that gi ~ 1 /S and (pi ~ VS for <5 — > 0). 
We can, furthermore, analytically integrate Eq. ED by 
neglecting the effect of fi, Vi and gi as compared to Eg 
(and by using the fact that noise amplitudes \i on adja¬ 
cent grid points are uncorrelated, in striking contrast to 
the coherent amplitudes ipi)- This yields 

OO 

Xi(t)= Y, MO)i l - l 'Ji-i'(Est/n)e- iE +/ h (33) 

l '=—oo 

where Ji denotes the Bessel function of the first kind of 
the order l. We thereby obtain the temporal correlation 
function of the quantum noise as 

xMxi'W) = {E s {f - t)/h) , 

(34) 

which reproduces Eq. ED as special case for t = t'. 

The third line of Eq. ED can then be rewritten as 
£z (t)ipi(t) +7with the effective (real and com¬ 
plex) noise terms 

m = gi(2\ X i(t)\ 2 - 1), (35) 

m(t) = gixf(t ) (36) 



Figure 1: Density and current profiles of the scattering state 
in the near-resonant transport configuration for two different 
grid spacings S = 0.5<5o (thick red lines) and 6 = 0.25o (thin 
blue lines). The upper panel shows the external potential 
V(x) (solid black line) as well as the interaction strength g(x ) 
(violet dashed line) according to Eqs. (1391 and (1401) . respec¬ 
tively. The green dotted line marks the level of the chemical 
potential. The middle panel shows the total densities (upper 
solid lines) and the incoherent parts of the density according 
to Eq. (1291) (lower dashed lines), while the lower panel shows 
the total currents, respectively for 8 = 0.5<5o (thick red lines) 
and 5 = 0.2<5o (thin blue lines). These tW calculations were 
done using 10® random realizations of the vacuum fluctua¬ 
tions. 


that satisfy £i(t) = rji(t) = r]i(t)r)i<(t ,s ) = 0 as well as 

&(*)&'(*') = 9i9i'Ji'-i (Es(t' — t)/H) , ( 37 ) 

»£ (*)•»'(*) = 

Without entering into too much technical detail, we can 
observe that we are dealing here with continuous-time 
stochastic processes with zero mean value whose ampli¬ 
tudes scale as gi/h and which vary on a characteristic 
time scale H/Eg. This yields effective diffusion constants 
that scale as {gi/h) 2 (h/E$) = mg 2 (l5)/h 3 and are there¬ 
fore finite in the continuous limit <5—^0. 

Note that the replacement \tpi\ 2 —> — 1 in 

the translation of the nonlinear term from the Gross- 
Pitaevskii to the tW equation is crucially important for 
the above argument to hold. In colloquial terms, the 
effective potential V) eff = —gi that is artificially intro¬ 
duced by this replacement exactly compensates for the 
enhanced repulsive interaction of the nonlinear wave with 
the vacuum fluctuations at positions where gi is particu¬ 
larly large. 

B. Numerical results 

Let us now verify these findings for the specific case of 
the transmission of a Bose-Einstein condensate through a 
symmetric double barrier potential. The latter is formed 
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Figure 2: Same as Fig. |T] for the non-resonant transport con¬ 
figuration. Total and incoherent densities (middle panel) as 
well as total currents (lower panel) were computed for the 
grid spacings <5 = 0.2739<5o (thick red lines) and 8 = 0.1095<5o 
(thin blue lines), using again 10 6 random realizations of the 
vacuum fluctuations. 


by a sequence of two Gaussian barriers, i.e. 


V(x) = Vo i) 2 /(2<r 2 ) + e -(x-x 2 )V(2-T 2 )j ) ( 39) 

which are placed at the positions X\ and X2 and which 
have the same width a. The spatial dependence of the 
one-dimensional interaction strength is chosen as 


g(x) = 9o 



X — X\ 
G 


— tanh 


X — X 2 
G 



(40) 


For the sake of simplicity, we consider the coupling 
strength to the trap to be strongly localized in space, 
such that we can write y/N k(x) = SqS(x). This source 
emits a matter-wave beam whose chemical potential /i is 
tuned such that it is of the same order as the potential 
barrier height. In the absence of the scattering poten¬ 
tial and the interaction, such a source would yield the 
homogeneous densities and currents 

Po = \S 0 \ 2 m/{2h 2 p) , (41) 

jo = \J2p/mp 0 ■ (42) 


We shall, in the following, consider two distinct trans¬ 
port scenarios: a near-resonant one with well-separated 
potential barriers, leading to nearly perfect transmission 
of the matter-wave beam, and a non-resonant one with 
overlapping barriers, leading to almost total reflection of 
the atoms. The former (near-resonant) configuration is 
obtained by choosing the parameters x\ = X\ = 25<5o, 
X 2 = X 2 = 75<5o, a = a = 56 q, go = 0.034Vb<5o, 
So = O.SVov^ and p = 1.15Vo- Here we introduce 
by 5 0 = h/y/mV o the characteristic length scale that 
corresponds to the barrier height Vo- For the latter 
(non-resonant) configuration, we choose X\ = 11.50^o, 
X 2 = 15.88i5o, er = 1.095<5o, x\ = 3.834<5o, i '2 = 24.1(5q, 


a = 1.643<5 0 , g 0 = 0.018261^0, S 0 = 1.974y 0 v^o, and 
p = 0.727Vb- 

The practical numerical implementation of the tW 
method in this waveguide scatteri ng p roblem is done in 
the same way as described in Ref. |24[. In particular, we 
restrict the numerical representation of the waveguide to 
a finite spatial region (0 < x < 100do in the near-resonant 
and 0 < x < 27.39(5o in the non-resonant case) in which 
the scattering potential and the interaction strength can 
be considered to be nonzero. Within the framework of 
the tW approach, the coupling to the two noninteracting 
“leads” on the left- and right-hand side gives then rise 
to a time-dependent quantum noise that enters into the 
scattering region |24[. The decay of atoms to the leads, 
on the other hand, is modeled through smooth exterior 
complex scaling |34| . 

In the numerical practice, the source amplitude is 
slowly increased with time from zero to its maximal value 
So- Figure Q] shows the stationary density and current 
profiles that are numerically obtained after this ramp¬ 
ing process in the near-resonant transport configuration, 
for two different choices for the numerical grid spacing 5. 
We clearly see a nearly symmetric density profile, which 
is a signature of near-resonant transmission. The den¬ 
sity maxima at the positions of the potential barriers as 
well as the slight enhancement of the density in the inter¬ 
acting region in between the barriers are explained by an 
effective decrease of the speed of the atoms due to energy 
conservation, which renders the atoms more likely to be 
detected there. As is indeed expected to occur in the 
presence of quasi-stationary scattering, the current pro¬ 
file, on the other hand, is fairly homogeneous within and 
outside the interacting region, apart from statistical fluc¬ 
tuations that arise from a finite Monte-Carlo sampling in 
the framework of the tW method. 

Figure^] shows the stationary density and current pro¬ 
files in the non-resonant transport configuration, again 
for two different choices for the numerical grid spacing 
8. While the chosen chemical potential p = 0.727Vo lies 
rather close to a single-particle resonance of the double 
barrier potential, the presence of the atom-atom interac¬ 
tion gives rise to a significant shift of the effective reso¬ 
nance level to higher values of the chemical potential and 
thereby induces blocking of resonant transmission [l5j . 
Again, a fairly homogeneous current profile is obtained 
within and outside the interacting region. 

Figure [3] shows how the total and incoherent transmis¬ 
sions across the double barrier potential scale with the 
grid spacing S. We calculated for this purpose the average 
total current j and the incoherent current j lncoh accord¬ 
ing to Eqs. m and ®, respectively, and performed 
an additional average over the numerical grid points of 
the system, in order to reduce the size of the statisti¬ 
cal fluctuations. Clearly, we see that both the total and 
the incoherent transmission tend to finite values in the 
limit of vanishing grid spacing <5 —> 0, both for the near¬ 
resonant and the non-resonant transport configuration. 
This numerical observation confirms that the tW method 
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Figure 3: Total (blue circles) and incoherent transmissions 
(red diamonds) across the double barrier potential (1391) . ob¬ 
tained with the tW method for different values of the numer¬ 
ical grid spacing 5 in the near-resonant (upper and middle 
panels) and the non-resonant transport configuration (lower 
panel) using the same parameters of the scattering system as 
for Figs, m and 0 respectively. The green solid lines show 
the results of a purely classical Gross-Pitaevskii calculation 
according to Ref. [1E|]. The total and incoherent currents 
are calculated according to Eqs. (12611 and (|30l>. respectively, 
where we perform an additional average over the numerical 
grid points of the system (which then defines the depicted 
error bars). In both transport configurations, the total and 
the incoherent transmission tend to finite values in the limit 
of vanishing grid spacing <5 —> 0 [35| . 


is expected to yield consistent results in the continuous 
limit. 


ing perspectives for the applicability of the tW method 
in the context of bosonic waveguide scattering processes. 
Indeed, judging from Ref. we expect that the tW 
method will yield reliable predictions in the presence of 
weak interaction strengths, similarly as for the Hartree- 
Fock-Bogoliubov method [19} or the Bogoliubov back- 
reaction method [36l . [37]. It certainly breaks down in 
the presence of very strong interactions where genuinely 
quantum many-body approaches based on the density 
matrix renormalization group (DMRG) and the matrix- 
product state (MPS) method [38-42] as well as on the 
Gutzwiller ansatz (see, e.g., Ref. |43[) are more appropri¬ 
ate to describe the dynamics of the system. 

We furthermore presented a semiclassical derivation of 
the tW approach in the framework of the van Vleck- 
Gutzwiller theory, which essentially relates this approach 
to the diagonal semiclassical approximation. This latter 
framework opens possibilities for a wider application of 
the tW method in the context of moderately interacting 
systems that exhibit chaotic classical dynamics. Indeed, 
while the conventional implementation of the tW method 
is not expected to yield reliable predictions for long evo¬ 
lution times that exceed the Ehrenfest time of such a 
chaotic system, average transport observables, which are 
obtained, e.g., in the presence of disorder, are never¬ 
theless expected to be correctly reproduced by this ap¬ 
proach provided we can safely assume classical ergodicity 
(and account for systematic quantum interference effects 
such as coherent backscattering in Fock space Ell). This 
new semiclassical perspective of the tW approach will 
also form a useful basis for the development of a truly 
semiclassical van Vleck-Gutzwiller (or Herman-Kluk) ap¬ 
proach that is able to capture interference effects beyond 
the diagonal approximation [2?],13: SI] ■ 


V. CONCLUSION 

In summary, we investigated the feasibility of the tW 
approach for the description of quasi-stationary scatter¬ 
ing processes with interacting Bose-Einstein condensates 
in one-dimensional waveguides. While a discretization of 
position space is most conveniently introduced in order 
to apply the method in practice, we showed that consis¬ 
tent results are to be obtained in the continuous limit 
of vanishing grid spacing. This finding provides promis- 
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